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A recently advanced argument against the atmospheric greenhouse effect is refuted. A planet with- 
out an infrared absorbing atmosphere is mathematically constrained to have an average temperature 
less than or equal to the effective radiating temperature. Observed parameters for Earth prove that 
without infrared absorption by the atmosphere, the average temperature of Earth's surface would 
be at least 33 K lower than what is observed. 

PACS numbers: 92.60.Vb,05.90.+m 

I. INTRODUCTION 

The results presented here are not new. However the form of presentation is designed to clearly and accurately 
respond to recent claims^ that a physics-based analysis can "falsify" the atmospheric greenhouse effect. In fact, the 
standard presentation in climatology textbooks^ is accurate in all material respects. The following explores in more 
detail certain points that seem to have been cause for confusion. 

First presented are the definitions of basic terms and the relevant equations for the flow of energy. The situation for 
a planet with no infrared-absorbing atmosphere is then examined, and a constraint on average temperature is proved. 

Several specific models of planets with no infrared-absorbing atmospehere are then solved, including one presented 
by Gerlich and TscheuschnerA, and it is verified that all satisfy this constraint. 

A simple infrared-absorbing atmospheric layer is added to these models, and it is proved that the temperature 
constraint is easily violated, as is shown by the observational data for Earth. 

II. DEFINITIONS AND BASIC EQUATIONS 

Define the incoming irradiance S as the energy per unit area and per unit time arriving at a planet from a stellar 
source. The actual radiation field is characterized by a spectrum of wavelengths and (depending on the size of the 
star(s) and distance to the planet) a small spread in directions. S is an integral over all wavelengths and propagation 
directions of the radiant specific energy at the distance of the planet from the star. As the planet moves through its 
yearly orbit, the value of S will vary, so it is strictly a function S(t) of time. 

The energy per unit time arriving at the planet is the product of S(t) with the area the planet subtends in the 
plane perpendicular to the radiant propagation direction. For a spherical planet of radius r, this area is simply irr 2 . 
So total arriving energy (energy per unit time, or power) from space is 

E in (t) = nr 2 S(t) (1) 

For a given location x on the planetary surface, the normal to the plane of the local surface makes an angle 6*(x, t) 
with the radiation propagation direction of a given stellar source. Only one side of the planet will be lit at any given 
time from that source; this can be generally indicated by those angles 9 from to ir/2 radians. Angles from ir/2 to ir 
would be unlit. So local incoming irradiance would be: 

s(x,i) = cos(0(x, t))S(i) (2) 

Integrating this on a spherical planet (including only the lit side) gives the irr 2 factor in equation 1. 

Define the albedo a of the planet as the fraction of incoming irradiance that is reflected, a is also a local property 
(for Earth much is reflected by clouds and ice), so the locally reflected energy is a(x, i)s(x, t). Integration across the 
lit side of the planet gives a well defined reflected energy: 

Ereflected{t) = S(t) J <x(x, t) COs(0(x, t))dx (3) 

An effective albedo a e ff can then be defined by the ratio of reflected to incoming energy across the planet as a 
whole: 
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E r eflected{t) = Cl e ff{t)E m (t) = TTr 2 a e ffS{t) (4) 

The difference between incoming and reflected energy is what the planet absorbs (again, per unit time): 

E a bsorbed{t) = ITT 2 (I - a e ff(t))S(t) (5) 

The fundamental characteristic of a planet in space is that of no material interchanges with its surroundings. The 
only substantive way energy can come in is through electromagnetic radiation, and the only way energy can leave 
is similarly through the planet's own electromagnetic emissions. There is a very small correction from gravitational 
tidal forces, and a planet also receives a small net energy input from internal radioactive decay, but for planets like 
the Earth these are thousands of times smaller than the stellar input. 

A planet with no incoming absorbed energy would reach thermodynamic equilibrium with the cosmic microwave 
background, with a uniform temperature of about 2 K. Absorption of incoming stellar irradiance results in heating 
until a steady state with equal incoming and outgoing energy (measured outside the atmosphere, and averaged over 
one planetary revolution, or whatever the most important variation in time) is reached. Define T(x, t) as the local 
surface temperature of the planet, and e(x, t) as the local emissivity. Thermal radiation from the surface is then given 
by the Stefan-Boltzmann law^: 

E emi tted{t) = a J e(x, t)T(x, tfd* (6) 

Similar to the effective albedo, an effective emissivity and effective radiative temperature can be defined as averages 
over the planetary surface: 

T e//(*) 4 = ^/nx,i) 4 dx (7) 

and 

£e/j( ^ ^f / e(X ' t)T(X ' f)4& (8) 
Total radiated thermal energy from the surface can then be written in terms of the effective temperature and albedo: 

E emltt ed{t) = 4irr 2 ae eff {t)T eff {t) 4 (9) 

For a planet with no atmosphere, or with an atmosphere that doesn't absorb electromagnetic radiation to any 
significant degree, all this surface-emitted thermal radiation escapes directly into space, just as all the absorbed 
stellar radiation reaches the ground. So the net rate of change in energy of the planet at time t is: 

E p lanet{t) = E a bsorbed{t) - E em itted{t) = TTT 2 (I - a e ff(t))S(t) - AtTT 2 <je e f f [t)T e f / (t) 4 (10) 

We will look at relevant constraints associated with an absorbing atmosphere via a simple model later in the 
discussion. 

The orbital processes for a planet (and any internal variability in the star) determine the variation in S(t); that 
combined with rotation and internal dynamics gives variability in a e ff(t), e e ft(t), and T e ff(t). As a result the planet 
may experience natural periods of warming or cooling as E p i anet (t) goes positive or negative, respectively. On average, 
however, over time, this rate of energy change should come very close to zero as long as all the input parameters are 
reasonable stable over the long term. If it didn't average to zero for a long period of time, the energy of the planet 
would cumulatively build or decline. 

In addition to the effective temperature obtained from averaging temperature to the fourth power, relevant for the 
thermal radiation problem, we should also look at the more natural average temperature for the planetary surface: 
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T ave (t) = — j T(x,t)dx (11) 

As Gerlich and Tscheuschner note 1 - in their section 3.7, thanks to Holder's inequality, this average temperature 
T aV e(t) is always less than or equal to the effective thermal radiation temperature T e ff(t), so T^ ve is less than or equal 
to T 4 ^ j , and rearranging Eq. [10] gives the following constraint on average temperature. 

T ave (t) 4 < ^tt((1 - aeff{t))S(t)/4 - E planet {t)/±izr 2 ) (12) 



III. SOME EXAMPLES 



A. Model 1: Nonrotating planet 



First let's look at the simple model planet solved by Gerlich and Tscheuschner (section 3.7.4) This is a non-rotating 
planet (or a planet with a rotation axis parallel to the incoming radiation) with no internal heat transport in constant 
local radiative equilibrium so that E p i anet is always zero. The non-rotation removes all time-dependences. Emissivity 
is assumed to be 1 everywhere; likewise S and a are uniform. Also the microwave background is ignored so the un-lit 
side of the planet is always at absolute zero temperature. From Eq. [2] for the local irradiance s(x) we quickly obtain 
the local temperature for the lit side of the planet: 



T roode n(x) - {(1 - a) cos(#(x))S/a} 1/4 (13) 
The average temperature is obtained by integrating over the sphere: 

T „ e =i((l-a)5/a) 1 / 4 r /2 cos(0) 1 / 4 2 7 rsin(0)^ = ^((l- a )5/a) 1 / 4 (14) 
4?r Jo 5 

The effective temperature similarly is given by 

i _ , r /2 i 



T eff = ^((1 - o)S/a) J co S (0)27r S in(0)d0 = -((1 - a)S/a) (15) 
which is what it has to be to ensure (Eq. I10p that E p i anet is zero. 

So in this case the ratio T ave /T e ff is 2V2/5 or about 0.566, and the planet's average temperature is indeed well 
below the effective temperature in this simple model. Plugging in numbers appropriate for Earth, T e tf comes to 255 
K and T ave would be 144 K, for this non-rotating atmosphere-free version of the planet. 



B. Model 2: Simple rotating planet 



To this simple model let us now add rotation, including a local heat capacity effect that accounts for some heat 
transport vertically, while still leaving out any transport of heat horizontally from one location to another. Assume 
the radiation direction is in the plane of rotation. Define the rotation period D (a day for the planet) and a thermal 
inertia coefficient^ c with units of J/Km 2 . On a real planet, c depends on temperature and on D (a time- or frequency- 
dependence); physically it represents the product of the volumetric heat capacity and the depth or height to which 
the incident heat energy is circulated or conducted during a daily thermal cycle, c then determines the local rate of 
change of temperature based on the local version of the net energy equation: 



cT(x,t) — E a bsorbed(x, ^emitted 

(x,t) = {1 ~ a)S cos(e(x,t))Q(cos{e{x,t))) - aT(x,t) 4 (16) 

Represent x by angular coordinates <fi for longitude and £ for latitude. Thanks to the rotation, the position of x 
relative to the incoming sunlight changes as if <fi were steadily incremented at a rate 2nt/D. The sun angle 9 then is 
found from: 
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cos(0) = cos(cj) + 2nt/D) cosf (17) 

The driving forces in the equation repeat with period D, so under steady state conditions the solution(s) of Eq. [TB] 
should also repeat with that period. At any point 4> 1 £ on the surface this solution would follow exactly the same curve 
of temperature as a function of time for every longitude 4>, at the given latitude £. In the following we replace 2tt / ' D 
with the symbol u and (1 — a)S cost; with A, and for simplicity set tj> = —ir/2 (for all other points on the surface the 
solution is just shifted forward or back a bit in time). Then the step function 0(cos(0 4- urf)) becomes equivalent to 
a square wave W(u)t) which is 1 for ujt between and w, for u>t between tt and 2ir, and repeating periodically after 
that. 

So Eq. [16] is reduced to: 



cT = A sin(wt) W(ut) ~ aT 4 (18) 

Any non-transient solution for T(t) will be periodic in time so that T(D) = T(0). Integrating Eq. [T8l over a 
planetary day (t = to t = D) gives: 



2A f D A 

= a I T 4 dt (19) 



Define an effective radiative temperature T e ff(£) for latitude £ based on the average fourth power (whether averaged 
over time or over longitudes is the same): 



T «ff(0 A = ^l T ' dt ( 2 °) 

then rearranging Eq. 1191 and substituting in the definitions of to and A gives: 

Te//(e)4 = d-^cosg (2i) 
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Note that the peak T e jf for £ = (the equator) is a factor of l/n 1 / 4 sa 0.75 times the peak temperature on the 
non-rotating planet (from the point directly under the sun). 

Once again we can check that the rate of change in net energy for the planet as a whole (Eq. [TO"]) comes to zero by 
finding the effective radiative temperature for the entire surface: 



1 r /2 (l-a)Scos£ „ ,„ (l-a)S r /2 2> „ (1 - a)S 

T -ff = T- / ' 27rcos ^ = V ^ / cos 2 ^ = V ; 22 

4vr./_ 7r/2 7TO- 2trt J_ 7r/2 4cr 



as it has to be, the same as for the non-rotating case in Eq. [TBI 

While we won't find a full analytic form for the temperature as a function of latitude and time in this model, we 
can learn a bit by examining the time dependence of temperature in Eq. [15] more closely. First, define x = cot and 
y(x) = T(£, t)/T e ff(£) for a given latitude £, with T e //(£) determined by Eq. [21] Eq. [T5]then reduces to: 



dy (1 — a)Scos£ , , . , ,. , crT e fr(£) 3 A , . . . . , , 1 d . 

/ = rJ lt , Mx)W(x) e -U^y 4 = X(sm(x)W(x) - -y 4 ) (23) 

where we define the dimensionless parameter A = (1 — a)S cos^/c toT e ff(Jf) Physically this roughly represents the 
ratio of the quantity of incoming energy absorbed in a day to the total heat content of the surface (to a relevant 
depth) at the effective radiative temperature. If A is small (heat capacity or rotation frequency high or latitude close 
to the poles), heating or cooling will occur only slowly, and the temperature will stay close to T e ff throughout the 
day (y will be close to 1). If A is large (heat capacity or rotation frequency low, latitude closer to the equator) then 
heating and cooling are rapid, and the temperature variation is more significant. 

We can find an analytic solution for the night side of the planet, where the square wave W(x) = (x between 
(2n — 1)tt and 2nir for integer n). Eq. [53] loses all dependence on x and is easily integrated: 



Temperature on a rotating planet 




FIG. 1: Temperature relative to the effective radiating temperature (T e //) for the simple rotating planet model for various 
values of the thermal response parameter A. The plots are of temperature against time where is sunrise, ir is sunset, and 2n 
is sunrise again. 
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y( x ) = ( 



3A(x — a) 



a/3 



(24) 



where a is a constant determined by the initial condition y(7r) (the value of the temperature when night begins for 
X = 7r): 



a = 7r(l 
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3Ay(^)3 



y(27r)/y(7r) = 1/(1 + 3AJ/W 3 ) 1 / 3 



(25) 



which gives us the night-time temperature drop. When A is small and y is not too large to start with, the change is 
small - a fractional decline of roughly Ay(7r) 3 . For large values of A, the night-time temperature drop is limited by this 
slow inverse 1/3 power. This makes sense as the rate of temperature decrease must decline sharply as temperature 
gets lower and the T 4 radiative term drops. 

For the daytime, if the temperature starts out low with y « 1, then the y 4 term is negligible, at least at first, and 
Eq. [221 can be integrated easily enough: 



dy 
dx 



A sin ; 



y(x) 6 — A cos x 



(26) 



where the inegration constant b again is set by the initial condition b — y(0) + A, which y(n) — y(0) + 2A, i.e. the 
temperature increments by 2 A on a sinusoidal curve during the day. Of course when A is large or y(0) starts close to 
one, this approximation breaks down. 
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FIG. 2: The average, minimum, and maximum values of the relative temperature y from numerically integrating the equations 
for different values of A. The numerically computed average value of y 4 is also shown; this should always be exactly 1. Lines 
are shown indicating the approximate values of A corresponding to the equators of the terrestrial planets, which gives a picture 
of their most extreme temperature profiles if they had infrared-transparent atmospheres. 

Another approximation is to assume variations in y are small and that we can linearize about some chosen value 
yo- This gives: 



dy , . 
— — w A sin x 

ax 



^-M(y-y ) 



(27) 



which as a linear ordinary differential equation yields a solution: 



A 



1 + (4Ay » ; 



■ sin x — cos x 



(28) 



(3 here is another constant of integration to be determined by an appropriate initial condition. The linearization 
fails once y deviates significantly from yo, but the result should be generally valid if A is small, and can be used to 
generate step-wise solutions for the daytime temperatures under any value of A; numerical intergration of the basic 
equation of course can do the same. 

Numerical computed solutions for various values of A are shown in Fig. [TJ Fig. [5] shows the trends for average y 
(T/T e ff), average fourth power of y, and the minimum and maximum y values as A increases. As expected, the average 
fourth power is fixed at 1, while the average y decreases as A increases, eventually approaching the non-rotating value 
of 2v2/5 as A — * oo. Also note the maximum temperature quickly approaches the nonrotating value of y = ypn for 
large A. The minimum temperature drops slowly as y ~ 1/(3A) 1 / 3 . 

Approximate formulas for the average temperature are, for large A: 



2V2/5 + 0.392A 



-0.279. 



A 



(29) 



and for small A: 



y me ~l-0.196A 2 ;A^0 



(30) 
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Note that from Eq. [21] the temperature scale varies with latitude as cos(^) 1 / 4 , while the value of A varies as 
cos(£) 3 / 4 . So finally, integrating over the whole planet we have an average temperature value of: 



T aw = C- ^)V4/ C os(£) 1/4 2w(A(£))cos(£R (31) 

™ Jo 

For small A we could substitute the expression from Eq. 1301 in any case, we know y ave < 1 for all latitudes, so we 
have an upper bound on the average temperature of the entire rotating planet T ave by substituting in the numerical 
value for the cos(£) 5 / 4 integral: 



T ave < 0.69921(^ ^)!/ 4 (32) 

a 

and note that this bound is a little over 1% less than T e ff for the entire planet (Eq. |2"2")) which has a constant 
(1/4) 1 / 4 = 0.7071... instead of 0.69921 in the same expression. 

So no matter the rotation rate, no matter the surface heat capacity, the average temperature of the planet in this 
rotating example, with only radiative energy flows and no absorbing layer in the atmosphere, is always less than the 
effective radiating temperature. For very slow rotation or low heat capacity it can be significantly less; for parameters 
in the other direction it can come as close as 1% (i.e. up to 252 K on a planet like Earth). 



C. Model 3: Rotating planet with varying albedo 



While the variability in infrared emissivity is relatively small across the surface of a realistic planet, the albedo can 
be significantly different from place to place. One of these involves taking into account the effect of ice, by which 
the high latitudes reflect more incoming radiation back into space than equatorial latitudes do. What effect does this 
have on effective radiating temperature and total temperature? 

We can model this by a slight change in model 2, by making the value of a dependent on latitude £. For example let 
a = sin 2 (£) so it is zero at the equator, and approaches 1 at the poles. This changes nothing in most of the analysis 
of the preceding section, until we integrate over latitudes. For Eq. [22] we now have: 

T* ff = — / ^ '- i.27rcos^=— / cos 4 ^£= — 33 

This just means that our sin 2 albedo has the same effect on the total radiation absorbed by the planet as would a 
uniform albedo value of 1/4. However, it redistributes that energy, putting more near the equator and less near the 
poles. The effect on average temperature across the planet, for this modified version of Eq. [3T1 is: 

T ave = ( — ) 4 / 4 cos(0 3/4 y a „ e (A(0) cos(£R (34) 

™ Jo 

which then, putting in the numerical value for the cos 7 / 4 integral, gives the inequality 



TABLE I: Relevant parameters for the planets. See |http://nssdcTg sfc. nasa.gov/planetary/factsheet/ A for Eq. [23] (at the 
equator, £ = 0) estimated from thermal inertia, solar day, and the other parameters. It is particularly small for Earth thanks 
to rapid rotation and the high heat capacity of water covering most of the surface. 
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T ave < 0.6206(-) 1 / 4 (35) 
a 

which is about 5% below the effective temperature (the numerical coefficient is (3/16) 1 / 4 or 0.6580). 

IV. INFRARED ABSORPTION IN THE ATMOSPHERE 

The examples of these simple models show that vertical energy transport for a planet with a transparent atmosphere 
only smooths out the daily temperature curve, without being able to bring the surface temperature higher than the 
effective radiative temperature. The same is true if we were to add in more realistic horizontal energy transport from 
larger-scale atmospheric and oceanic circulation - of course getting much more realistic means entering the realm of 
more full-scale general circulation models^, which we have no intention of doing here. 

On a planet with significant internal energy sources the effective temperature for radiative balance could be exceeded 
even with a transparent atmosphere. For example a planet still losing its initial heat of formation, or a planet remote 
from its sun with a high enough radioactive content, or on a planet or moon with very large tidal forces, you will 
have a net outward flow of energy to space, and may well have an average temperature above the limit. But for the 
terrestrial worlds of our solar system, these internal sources of heat are thousands of times too small to have any 
noticeable effect on surface temperature. 

And yet the observed average surface temperature on Earth and Venus significantly exceeds the effective radiative 
temperature set by the incoming solar radiation. This is not observed for the Moon or Mars. What makes Venus and 
Earth so different? 

Net energy flux is determined by the radiation that gets into space, not what leaves the surface. The only way for a 
planet to be radiatively warmer than the incoming sunlight allows is for some of that thermal radiation to be blocked 
from leaving. That means some layer above the surface must be absorbing or reflecting a significant fraction of the 
outgoing infrared radiation. I.e. the atmosphere must not be transparent to infrared. 

So, let's add to our rotating planet model a simple model of this blocking effect: a fraction / (between and 1) of 
the outgoing radiation E em itted from the surface is absorbed by a thin layer of the atmosphere. This layer will have 
its own temperature but for simplicity we make the assumption that the heat capacity of the atmospheric layer is low 
so that it remains essentially radiatively balanced through the day, and the specific temperature becomes irrelevant. 
That means that this atmospheric layer continuously emits an amount / • E em i t t e d equal to what it absorbs from the 
ground. 

Since thermal re-emission is randomly directed, half the radiation from this atmospheric layer will go up, and half 
down. Assuming the surface is fully absorbing and the rest of the atmosphere is transparent, total outgoing radiation 
from the planet (above the atmospheric layer) is then: 

Eout = (1 — f)E em itted + —j 'E em itted = (1 — f /fyE emitted (36) 

while absorbed radiation on the surface is now: 

E a bsorbed{t) = 7Tr 2 (l - a e ff(t))S(t) + -j 'E emiU ed (37) 

Incoming solar radiation still drives everything - if the solar constant S drops, then so does everything else. But 
the effect of the absorbing layer is to reduce the final outgoing energy for a given temperature, so the planet heats up 
until things are back in balance again. 

Generalizing Eq. [TO] we have net energy change (which can be calculated either at the surface or above the absorbing 
layer of the atmosphere): 

E surface (t) = E absorbed (t) - E emitted (t) = - 2 (1 - a eff (t))S(t) + \fE emMed {t) - ^ae e}f (t)T ef} {tf (38) 

= vrr 2 (l - a eff (t))S(t) - Avr 2 oe eff (t)(l - //2)T e// (t) 4 (39) 

We then end up with essentially the same equations as in the previous section, for example Eq. [TO] is the same, 
except that effectively the solar input S and thermal inertia c in those equations are increased by the factor 1/ (1 — //2). 
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That means the surface effective radiative temperature T e tf in those equations is increased by a factor (1/(1 — 
//2)) 1 / 4 , or as much as 2 1 / 4 for a fully absorbing atmospheric layer. The parameter A is then reduced by that same 
ratio (the increases in S and c cancel out, leaving a 1/T e // term). So the temperature curve of this radiatively 
insulated planet is even more uniform than without the insulating layer. The average temperature can come within 
a few percent of this higher T e f f , or well above the limits for a planet with a transparent atmosphere. 

A more realistic atmosphere would be characterized by more than one absorbing layer (or a thick layer with a 
temperature differential and limited conductivity from bottom to top), which will further decrease outgoing thermal 
radiation and increase surface temperatures. Details of absorption in the real atmosphere also depend on pressure; 
nevertheless, the presence of any absorption at all is what qualitatively distinguishes a greenhouse-effect planet from 
one with a transparent atmosphere, and is what allows surface temperatures to climb above the effective radiative 
limit. 



V. CONCLUSION 

Gerlich and Tscheuschner— state, among more extravagant claims, that "Unfortunately, there is no source in the 
literature, where the greenhouse effect is introduced in harmony with the scientific standards of theoretical physics." 
The above analysis I believe completely establishes, within perfectly simple and appropriate theoretical physics con- 
structs, the main points. Namely that assuming "the atmosphere is transparent for visible light but opaque for 
infrared radiation" leads to "a warming of the Earth's surface" relative to firm limits established by basic physical 
principles of energy conservation, for the case of an atmosphere transparent to both visible and infrared. 

In particular, it has been shown that: 

1. An average surface temperature for a planet is perfectly well defined with or without rotation, and with or 
without infrared absorbing gases 

2. This average temperature is mathematically constrained to be less than the fourth root of the average fourth 
power of the temperature, and can in some circumstances (a planet with no or very slow rotation, and low 
surface thermal inertia) be much less 

3. For a planet with no infrared absorbing or reflecting layer above the surface (and no significant flux of internal 
energy), the fourth power of the surface temperature always eventually averages to a value determined by the 
incoming stellar energy flux and relevant reflectivity and emissivity parameters. 

4. The only way the fourth power of the surface temperature can exceed this limit is to be covered by an atmosphere 
that is at least partially opaque to infrared radiation. This is the atmospheric greenhouse effect. 

5. The measured average temperature of Earth's surface is 33 degrees C higher than the limit determined by items 
(2) and (3). Therefore, Earth is proved to have a greenhouse effect of at least 33K. 

The specific contributions of individual gases such as CO2 to Earth's greenhouse effect are covered well by the 
standard treatments of the subject^^. 



* E-mail: apsmith@aps.org 
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